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A mechanism explaining a strong enhancement of nonleptonic weak decays was 
\f\ ' suggested in 1975, later to be dubbed the penguin. This mechanism extends Wil- 

son's ideas about the operator product expansion at short distances and reveals an 
intricate interplay of subtle features of the theory such as heavy quark masses in 
Glashow-Iliopoulos-Maini cancellation, light quarks shaping the chiral properties 



of QCD, etc. The penguins have subsequently evolved to play a role in a variety 
of fields in present-day particle phenomenology. I will describe the history of this 
idea and review its subsequent development. The recent measurement of direct 
CP violation in K decays gives a new confirmation of the penguin mechanism. 



1 History of the idea 
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It was an exciting period, with Quantum Chromodynamics (QCD) emerging 
as the theory of strong interactions, when three of us - Valya Zakharov, Misha 
Shifman and I - started in 1973 to work on QCD effects in weak processes. 
qj ' The most dramatic signature of strong interactions in these processes is the so 

t£^ • called AI =1/2 rule in nonleptonic weak decays of strange particles. Let me 

remind you what this rule means by presenting the experimental value for the 
ratio of the widths of K$ — > tt + it~ and K + — ► 7r + 7r° decays 

X 

JH ' T(K S ^tt+tt-) 

The isotopic spin / of hadronic states is changed by 1/2 in the K$ — * tt + it~ 
weak transition and by 3/2 in K + — ► 7r + 7r°, so the AI = 1/2 dominance is 
evident. 

What does theory predict? The weak interaction has a current x current 
form. Based on this, Julian Schwinger suggested!!! to estimate nonleptonic 
amplitudes as a product of matrix elements of currents, i.e. as a product of 
semileptonic amplitudes. This approximation, which implies that the strong 
interaction does not affect the form of the weak nonleptonic interaction, gives 
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9/4 for the ratio (fil). Thus, the theory is off by a factor of two hundred! We 
see that strong interactions crucially affect nonleptonic weak transitions. 

The conceptual explanation was suggested] by Kenneth Wilson in the con- 
text of the Operator Product Expansion (OPE) which he introduced^ Assum- 
ing scaling for OPE coefficients at short distances, Wilson related the enhance- 
ment of the AI = 1/2 part of the interaction with its more singular behavior 
at short distances as compared with the AI = 3/2 part. In the pre-QCD era 
it was difficult to test this idea having no real theory of the strong interaction. 
With the advent of QCD all this changed. The phenomenon of asymptotic 
freedom gives full theoretical control of short distances. Note in passing that 
the American Physical Society also followed this development: the discoverers 
of asymptotic freedom, Gross, Politzer, and Wilczek, became recipients of the 
1986 Sakurai Prize. 

In QCD the notion of OPE in application to nonleptonic weak interac- 
tions can be quantified in more technical terms as one can calculate the ef- 
fective Hamiltonian for weak transitions at short distances. The weak inter- 
actions are carried by W bosons, so the characteristic distances are ~ 1/mw, 
with mw = 80 GeV. The QCD analysis at these distances in tha effective 
Hamiltonian was done in 1974 by Mary K. Gaillard with Ben Leeo, and by 
Guido Altarelli with Luciano MianiO. Asymptotic freedom at short distances 
means that the strong interaction effects have a logarithmic dependence on mo- 
mentum rather than power-like behavior, as was assumed in Wilson's original 
analysis. The theoretical parameter determining the effect is log(miy/AQCD) 
where Aqcd is a hadronic scale. These pioneering works brought both good 
and bad news. The good news was that, indeed, strong interactions at short 
distances logarithmically enhance the AI = 1/2 transitions and suppress the 
AI = 3/2 ones. The bad news was that quantitatively the effect fell short of 
an explanation of the ratio (1). 

Besides 1/mw an d 1/Aqcd there are scales provided by masses of heavy 
quarks t, b and c. In 1975 the object of our study was distances of order l/m c 
- the top and bottom quarks were not yet discovered. Introduction of top and 
bottom quarks practically does not affect the AS = 1 nonleptonic transitions. 
However, it is different with charm. At first sight, the c quark loops looked 
to be unimportant for nonleptonic decays of strange particles in view of the 
famous Glashow-Illiopoulos-Miani cancellationQ (GIM) with corresponding up 
quark loops. In 1975 the belief that this cancellation produced the suppression 
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factor 
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was universal, which is the reason why the effect of heavy quarks was over- 
looked. We found instead that: 

(i) The cancellation is distance dependent. Denoting r = 1/fi, we have 

for m c <C u < mw ; 

M 

o 
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log — f- , for fi < m c . (2) 

TVo suppression below m c \ 

(ii) Moreover, new operators appearing in the effective Hamiltonian at dis- 
tances larger than \/m c are qualitatively different - they contain right- 
handed light quark fields in contrast to the purely left-handed structures 
at distances much smaller than \/m c (see the next section for their ex- 
plicit form) . It was surprising that right-handed quarks become strongly 
involved in weak interactions in the Standard Model with its left-handed 
weak currents. The right-handed quarks are coupled via gluons which 
carry no isospin; for this reason new operators contribute to A/ = 1/2 
transitions only. 

(iii) For the mechanism we suggested it was crucial that the matrix elements 
of novel operators were much larger than those for purely left-handed 
operators. The enhancement appears via the ratio 

o 

2GeV, 



m u + m d 

which is large due to the small light quark masses. The small values of 
these masses was a new idea at the time, advocated in 1974 by Heiri 
Leutwyleid and Murray Gell-MannE The origin of this large scale is not 
clear to this day but it shows that in the world of light hadrons there is, 
besides the evident momentum scale Aqcd, some other scale, numerically 
much larger. 

Thus, the explanation of the AI = 1/2 enhancement comes as a nontriv- 
ial interplay of OPE, GIM cancellation, the heavy quark scale, and different 



intrinsic scales in light hadrons. I will discuss the construction in more de- 
tail below. However, I will first digress to explain what the mechanism we 
suggested has in common with penguins. 

We had a hard time communicating our idea to the world. Our first 
publication was a short letter published on July 20, 1975 in the Letters to 
the Journal of Theoretical and Experimental Physics. Although an English 
translation of JETP Letters was available in the West we sent a more detailed 
version to Nuclear Physics shortly after. What happened then was a long 
fight for publication; we answered numerous referee reports - our paper was 
considered by quite a number of experts. The main obstacle for referees was 
to overcome their conviction about the GIM suppression (mj: — rnu)/ m w an( ^ 
to realize that there is no such suppression at distances larger than l/m c . 
Probably our presentation was too concise for them to follow. 

Eventually the paper was published in the March 1977 issue of Nuclear 
Physics without any revision, but only after we appealed to David Gross who 
was then on the editorial board. The process took more than a year and 
a half! We were so exhausted by this fight that we decided to send our next 
publication containing a detailed theory of nonlcptonic decays to Soviet Physics 
JETP instead of an international journal. Lev Okun negotiated a special deal 
for us with the editor Evgenii Lifshitz to submit the paper of a size almost 
twice the existing limit in the journal - paper was in short supply as was 
almost everything in the USSR. We paid our price: the paper was published 
in 1977, but even years later many theorists referred to our preprints of the 
paper without mentioning the journal publication. 

On a personal note, let me mention that a significant part of the JETP pa- 
per was done over the phone line - Valya and Misha worked in ITEP, Moscow, 
while I was in the Budker Institute of Nuclear Physics, Novosibirsk. The phone 
connection was not very good and we paid terrific phone bills out of our own 
pockets. 

I think, that in recognition of our work in the world at large it was Mary 
K. Gaillard who first broke the ice - she mentioned the idea in one of her 
review talks. Moreover, she collaborated with John Ellis, Dimitri Nanopoulos, 
and Serge Rudaz in the workEj in which they applied a similar mechanism to 
B physics. It is in this work that the mechanism was christened the penguin. 

How come? Figure 1 shows the key Feynman diagram for the new opera- 
tors in the form we drew it in our original publications EL It does not look at 
all penguin-like, right? Now look how a similar diagram is drawn in the paper 
of the four authors mentioned above. 

You see that that some measures were taken to make the diagram reminiscent of 
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Figure 1: Appearance of new operators due to quark loops 
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Figure 2: Quark loops in B decays 



a penguin. Let me refer here to John Ellis' recollections|j on how it happened: 

" Mary K, Dimitri and I first got interested in what are now 
called penguin diagrams while we were studying CP violation in 
the Standard Model in 1976 ... The penguin name came in 1977, 
as follows. 

In the spring of 1977, Mike Chanowitz, Mary K and I wrote a 
paper on GUTs predicting the b quark mass before it was found. 
When it was found a few weeks later, Mary K, Dimitri, Serge Rudaz 
and I immediately started working on its phenomenology. That 
summer, there was a student at CERN, Melissa Franklin who is 
now an experimentalist at Harvard. One evening, she, I and Serge 
went to a pub, and she and I started a game of darts. We made a 
bet that if I lost I had to put the word penguin into my next paper. 
She actually left the darts game before the end, and was replaced 
by Serge, who beat me. Nevertheless, I felt obligated to carry out 
the conditions of the bet. 

E 



c John sent his recollections to Misha Shifman in 1995, who published them in the prefao 
to his book. 



For some time, it was not clear to me how to get the word into 
this b quark paper that we were writing at the time. Then, one 
evening, after working at CERN, I stopped on my way back to my 
apartment to visit some friends living in Meyrin, where I smoked 
some illegal substance. Later, when I got back to my apartment 
and continued working on our paper, I had a sudden flash that the 
famous diagrams look like penguins. So we put the name into our 
paper, and the rest, as they say, is history." 

I learned some extra details of the story from Serge Rudaz who is my 
Minnesota colleague now. He recollects that for him to beat John in darts was 
a miraculous event. John was a very strong player and had his own set of darts 
which he brought to the pub. 

2 Effective Hamiltonian 

Application of Wilson's OPE to nonleptonic AS = 1 decays means the con- 
struction of the effective Hamiltonian H as a sum over local operators Oi, 

F ei V) = \/2G F CVwX>(M)<W- (3) 
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Here V us , V uc i are elements of the Cabibbo-Kobayashi-Maskawa mixing matrix 
and /x denotes the so called normalization point, which is the inverse of the 
shortest distance for which the effective Hamiltonian is to be applied, Oi are 
gauge invariant local operators made out of quark and gluon fields, and a are 
OPE coefficients (c numbers). 

At distances larger than l/m c , i.e. at fi < m c , only light u,d,s quarks 
and gluons remain as building material for the Oi. Operators can be ordered 
according to their canonical dimension d, and the corresponding OPE coeffi- 
cients are proportional to (1 / 'mw) d ~ A ■ As a selection criterion for operators we 
used their transformation features in the limit of chiral SU(2)ixSU(2)fl sym- 
metry, picking up operators which are SU(2)/j singlets. Under this criterion, 
the operator of lowest dimension (d = 5) is of gluomagnetic type (magnetic 
penguins) , 

T = iSRa^XdL G« v , (4) 

where G" is the gluon field strength tensor, and t a (a = 1, ... 8) are 3x3 gen- 
erators of the color SU(3). The corresponding OPE coefficient is proportional 
to the strange quark mass m s . For this reason the magnetic penguins turn out 
to be not important in AS* = 1 transitions. They are important, however, for 
the b quark whose mass is large. We will return to this point later. 



The operator basis of SU(3)a, invariant operators of d = 6 consists of six 
four-fermion operators. The first four operators are constructed from left- 
handed quarks (and their antiparticles, which are right-handed), 

01 = SL^dLUL^UL- SLjpULULj^dL, (8f , Ai" = 1/2), 

02 = SLl^dh ULl^UL + SLJiiUL ULj^dL + 2SLliJ,dL dL^dt 

+2sL7nd L sli^sl, (8 d , AI = 1/2), 

03 = SLl^dh ULl^UL + SLlitUL ULj^dL + 2sL"/ t j,dL d^dh 

-ZsLl^dL sl^sl, (27, AI = 1/2), 

04 = SLl^dh ULl^UL + SLJiiUL ULj^dL 

-s L ^d L d L ^d L , (27, AI = 3/2). (5) 

Every quark field is a color triplet q l and summation over color indices is 
implied, <727 M <7i = (<&)* l^iliY- What is marked in the brackets are the SU(3) 
and isospin features of the operators. 

Two more four-fermion operators entering the set contain also right-handed 
quarks (in SU(3)_r singlet form), 

05 = s Llll t a d L (u Rl n a u R + d Rl H a d R + Sjan a s R ) , (8, AI = 1/2), 

06 = SLl»d L {uicfuR + dsrfdR + s R Ys R ) , (8, AI = 1/2). (6) 

Operators 05 and Oq are different by color flow only. 

The operator basis we have introduced is recognized now (some_doubts 
were expressed in the literature at the beginning). The standard seilid used 
presents some linear combinations of 0i_6- Actually, the set is five instead of 
six combinations, the completness of the basis was lost on the way, although 
it is not important within the Standard Model. 

2.1 Evolution 

The effective Hamiltonian (2) may remind the reader of the Fermi theory of 
beta decay with its numerous variants for four-fermion operators. While in 
many respects the analogy makes sense, the difference is that the standard 
model together with QCD allows us to fix all coefficients Cj. In the leading 
logarithmic approximation the evolution of the effective Hamiltonian at /i > m c 
was found in Refs.utl. Penguins do not appear in this range and the result for 



H (m c ) has a simple form: 
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where as{fj) is the running coupling 
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in the range with TV "active" flavors. The modification due to the b quark 
is not significant, a few percent numerically, and the t quark effects do not 
appear in the leading logarithmic approximation. 

Penguin operators 05,6 show up due to evolution at p, below m c , 
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where the anomalous dimension matrix p is 
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Numerical results for the OPE coefficients depend on the normalization 
point jjl, pushing /1 as low as possible maximizes the effect of the evolution. 
We chose the lowest jjl as the point where Qs(m) = 1- The values of ci.2,3.4 
are relatively stable, say, under variation of the c quark mass but the penguin 
coefficients 05,6 depend on it rather strongly. This is not surprising, of course, 
since penguins are generated in the interval of virtual momenta between jjl and 



m c . Numerically, even for /i as low as 200 MeV, the coefficients 05^ are rather 
small. Our 1975 estimates for C5 were in the interval 

c 5 =0.06 -=-0.14. (11) 

With the present day coupling as(mz) = 0.115, it would be about twice 
smaller. The values of the coefficients 05^ are small and unstable. We will 
return to the problem of OPE coefficients at the low normalization point in 
connection with the procedure of calculating matrix elements, essentially it is 
about matching for perturbative and nonperturbative effects. 

We also found the coefficient for the magnetic penguin operator (||) . Two- 
loop calculations were necessary for the purpose. Nowadays due to the efforts 
of few groups the next to leading approximation has been found for all types of 
weak processes, see, e.g., the reviewEJ. Although it is nice to have an accurate 
values of the OPE coefficients the main phenomenological effects come from 
matrix elements as we will see in the next Section. The theoretical uncertainty 
is much bigger there. 

3 Matrix elements and phenomenology of nonleptonic decays 

We used the naive quark model to find matrix elements of four-fermion oper- 
ators 0i_6 • The model implies a factorization for amplitudes of K mesons 
decays, in the hyperon decays it is also the case for all operators but 0\. It is 
clear that the factorization does not work at the range of /z where the strong 
coupling o.s{p) is small being in contradiction with a calculable evolution. For 
this reason, if the naive quark model works at all it is only at very low values 
of /1 where the evolution is complete, i.e. where as(fJ-) ~ 1. But then the 
theoretical accuracy of the perturbative OPE coefficients is not good because 
it is governed by the same as(fJ-)- 

So we employed the following strategy: let us use the naive quark model for 
hadrons together with a factorization somewhere at small \i but let us allow 
adjustment of the OPE coefficients from phenomenological fits. It does not 
involve many parameters. Predominantly three coefficients, c±, C4 and C5, plus 
a few nonfactorizable matrix elements in hyperon decays determine the bulk 
of nonleptonic amplitudes. New relations arising from such a fit are in good 
agreement with experimental data. 

3.1 AI = 3/2 transitions 

Let us start with the decay K + — > 7r + 7r°. It is a AI = 3/2 transition and its 
amplitude is determined by the matrix element of O4, 

M(K+ -> ^+7r°) = s/2G F V* B V u d ^ + A^Oi\K + ) . (12) 



In the valence quark model the factorization of this matrix element is visible 
from the Fcynman diagrams presented in Fig. 3. 
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Figure 3: Quark diagrams for K meson decays 

Consider, for instance, the diagram a, 

Mt^{AuLl^u L \Q)^ + \- SL ^d L \K+)^A{nL\l^ULy\0)^ + \(- SL ) 3 r(dLy\K + ) 
= T i {AuLl»u L \Q){K+\s Ll »d L \K+) , (13) 

where out of three terms entering the definition of O4, the first one factorizes, 
the second factorizes after Fierz transformation, and the third does not con- 
tribute. The matrix elements in Eq. (13) arc known from semilcptonic 7r — > [iv 
and K — ► nev transitions, 
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Accounting for all diagrams in Fig. 3 in a similar way we get 

M(K+ -» tt+it ) = zc 4 G F V: s V ud m\U ■ 
Comparing it with the experimental value 

|M(ir + ^7r + 7r°)| = 0.05 G F m\m^ 

I V / I CX p J\ " 



(14) 
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wc find 



c 4 w 0.25 . 



(17) 



what is about 1.6 times less that the theoretical estimate of C4. 

The consistency check comes from the AI — 3/2 hyperon decays. In Fig. 4 
quark diagrams for the A — ► pit~ decay are presented. The symmetry of wave 
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Figure 4: Quark diagrams for hyperon decays 



functions and operators under permutations of color indices (first discussed 
by Pati and Wooi-H) is important for the analysis, Namely, the operator 0\ is 
antisymmetric under permutation of quark color indices, operators 02,3,4 are 
symmetric, and operators 05,6 have no specific symmetry, in these quarks differ 
by their helicities. Baryon wave functions are antisymmetric in color, so only 
antisymmetric operators survive in diagrams a,b,c. For the symmetric A7 = 
3/2 operator O4 only the diagram c remains. For this diagram factorization 
takes place, 



(7r-p|0||A) = -^=(^ n|O 4 |A) = ^(7r-|d L7MUL |0)(p|i2L7^L|A) 
~ — y=Uq^u p ( 7 M + g 5.47^75 J u A . 



(18) 
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Using the value (IL7J) for the coefficient C4 we get predictions for the AI = 
3/2 hyperon decay amplitudes. They are collected in the Table 1, where the s 
and p wave amplitudes A and B are defined by 

M=-iG F mluf(A + Bj 5 )ui (19) 

We see that the predictions for s waves are in a reasonable agreement with 
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Table 1: AI = 3/2 amplitudes in hyperon decays 



the data, for p waves the experimental accuracy is too low for a conclusion. 



3.2 AI =1/2 transitions 

The above analysis of color symmetry and factorization in quark diagrams in 
Fig. 4 shows that the operators 0\ (diagrams a, b, and c) and 0\ (the diagram 
d) are dominant in hyperon decays. Matrix elements of 0\ are factorizable but 
that of 0[ are not. 

There exists a combination of amplitudes, 
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which does not contain 0[. Thus, the ratio of p and s waves for this combi- 
nation is predicted, 
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Experimentally this value is 33±10. The uncertainty can be reduced if the Lee- 
Sugawara relation 2SZ + A~t — VSEq = is used. Then the value is 27.4±2.5. 
Thus, we see an experimental confirmation of our description. 
The comparison of s wave A 
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with the corresponding experimental value —0.51 ± 0.10 fixes 
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(23) 



which is about four times larger than the theoretical estimate. 

The operators 05,6 are dominant in K decays. The naive quark model, 
together with a meson (m a fa 700 MeV) dominance in the s wave tttt channel, 
gives for the matrix element 
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The result for the amplitude of K$ — > 7r + 7r decay 

(tt+tt- \H cS \K s ) = (0.85 + 0.20) i G F m\m v , 



(24) 



(25) 



where the first number comes from ©5,6 and the second from Oi_4, matches 
the experimental value. 

4 Further developments 

4-1 Limit of large N co \ or , chiral loops 

As we see in the approach presented above the main uncertainty comes from the 
range of momenta of the order of hadronic scales. To make a consistent analysis 
in this range Bardeen, Buras, and Gerard suggested to use a description in 
terms of chiral meson dynamicsli3. It means that QCD is used to fix the 
effective Hamiltonian at /1 below charm threshold but large enough to have 
ag(/j) <C 1, particularly, they choose /j, ~ lGeV. Matrix elements of this 
effective Hamiltonian which accounts for momenta below 1 GeV are calculated 
with mesonic loops instead of quark loops. The use of chiral meson loops can 
be parametrically justified in the limit of large number of colors _/V c . 

Notice that penguin operators appear due to mixing with the left-handed 
operators Oi_4, this mixing is suppressed as 1/N C . In the dual mesonic picture 
it is hadronic vertices which bring in 1/N C . The smallness of 1/N C emphasizes 
once more the necessity of a large hadronic scale in the problem, it was m^jm q 
in the naive quark model. In chiral dynamics AI = 1/2 amplitudes also come 
out enhanced what confirm the penguin explanations of AI = 1/2. The chiral 
loops also improve predictions for AI = 3/2 transitions providing an additional 
suppression. 
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4-2 Penguin decays of B mesons 

We mentioned that "magnetic" penguins, i.e. the d — 5 operators of the kind 
given by Eq. (4) are of particular importance for B decays (see reviewed) . The 
Feynman diagram for b — > sg(j) transitions is presented in Fig. 5. In difference 

w 




g(i) 



Figure 5: Magnetic penguins for b — ► sg(j) transitions 

with strange particle decays there is no suppression due to a small quark mass 
(mb enters instead of m s ). From the diagram is seen that "magnetic" penguins 
are sensitive to heavy quarks in the loop and, for this reason, to possible 
deviations from the Standard Model prediction, 

Br(6 -> 87) = (3.5 ± 0.3) • 1(T 4 (26) 

The inclusive rate of B — > X s j decays was measured by CLEO and ALEPH 
collaborations, 

CLEO: Br(B -> X sl ) = (2.32 ± 0.57 ± 0.35) • 10" 4 
ALEPH: (3.11 ± 0.80 ± 0.72) • 10~ 4 

No deviation was observed. This CLEO experiment was recognized at this 
meeting, Ed Thorndike received the 1999 Panofsky Prize. 

4-3 Direct CP violation in K decays 

At this meeting Robert Tschirhart reported a_new measurement of direct CP 
violation in K decays by KTeV collaboration!!!] 

e i 

- = (2.8 ±0.41)- 10~ 3 , (27) 

e 

which is in agreement with the CERN NA31 experiment but different from the 
previous measurement in the Fermilab E731 experiment. 
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The direct CP violation is a crucial test of the Standard Model. Does the 
result fit Standard Model? In his presentation Tschirhart defined an answer to 
this question as debatable. Indeed, in the recent workEl authors claim that the 
Standard Model leads to the predictions substantially (few times) lower then 
the experimental result (27), and only for extreme values of input parameters 
the theory can be consistent with the data. 

My statement on the issue is that our theory of nonleptonic decays nat- 
urally predicts e'/e consistent with Eq. (27). Leaving detail for discussion 
elsewhereEj, let me make few remarks. 

Actually, relatively large paLues of e'/e within_such approach was obtained 
long ago by Gilman and WiseEEl, and by Voloshircil, but then the approach was 
unfairly abandoned. The argumentation was, see e.g. Ref.LS, that the operator 
05 enters with the small coefficient at /i = m c or /i = 1 GeV, so other operators 
are important. This criticism, however, is not relevant to the dominance of the 
operator O5 in the low normalization point where factorization takes place. 

In one loop approximation direct CP violation shows up as an imaginary 
part of the coefficient C5 

Im(V c * V cd ) a s (jn w ) rri^ 
Imc 5 (m c ) = -^- ^m^ « 0.12InuT cs ^) 

due to the diagram of Fig. 1 with the c quark in the loop (t quark contribution 
is small). In difference with CP even part, i.e. Recs, which comes from virtual 
momenta between /i and m c (thanks to the GIM cancellation), the CP odd 
part, Imcs, comes from a larger range between m c and mw- The corrections 
due to logarithmic evolution can be simply accounted for, they increase the 
value of Imc5(ro c ). Additionally Im 05(11) is increased by an evolution down 
to the normalization fi where a(ii) ~ 1 and factorization is applied. 

Accounting for the phenomenological value (23) of CP even part we find 
e'/e in ballpark of the data. Besides confirmation of the quark mixing nature 
of CP violation in the Standard Model this serves, somewhat surprisingly, as 
one more confirmation of our mechanism for AI = 1/2 enhancement. It is nice 
to get such a surprise on the eve of the Sakurai Prize. 

5 Conclusions 

Summarizing the development of twenty four years I cannot, unfortunately, say 
that theoretical understanding of AI — 1/2 is very much advanced. Progress of 
last years in technically difficult calculations of higher loops corrections to OPE 
coefficients is not crucial when the effect comes from large numbers in matrix 
elements. As we discussed above the enhancement of AI = 1/2 reflects the 
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existence of the large momentum scale in light hadrons. In case of glueballs the 
scale could be even larger, as it discussed in the next talk by Valya Zakharo-\£3. 
Moreover, the momentum scale in light hadrons related with this enhancement 
is so large that treatment of the c quark as heavy becomes questionable. In 
this sense the c quark is not heavy enough (let me remind, in passing, that a 
low upper limit for its mass was theoretically founder in the Standard Model 
before the experimental discovery of the c quark). 

Thus, we still have some distance to go so I finish by the sentence: 
Penguins spread out but have not landed yet. 
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